Fermi liquid theory for heavy fermion superconductors without inversion symmetry 
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We present the microscopic Fermi liquid theory for magnetic properties and transport phenomena 
in interacting electron systems without inversion symmetry both in the normal state and in the 
superconducting state. Our argument is mainly focused on the application to noncentrosymmetric 
heavy fermion superconductors. The transport coefficients for the anomalous Hall effect, the thermal 
anomalous Hall effect, the spin Hall effect, and magnetoelectric effects, of which the existence is a 
remarkable feature of parity violation, are obtained by taking into account electron correlation effects 
in a formally exact way. Moreover, we demonstrate that the temperature dependence of the spin 
susceptibility which consists of the Pauli term and van-Vleck-like term seriously depends on the 
details of the electronic structure. We give a possible explanation for the recent experimental result 
obtained by the NMR measurement of CePtsSi [Yogi et al.: J. Phys. Soc. Jpn.75 (2006) Of 3709], 
which indicates no change of the Knight shift below the superconducting transition temperature for 
any directions of a magnetic field. 



I. INTRODUCTION 

For a large class of superconductors, there exists in- 
version symmetry which allows the classification of the 
Cooper pair according to parity; i.e. a spin singlet 
or triplet pairing state realizes. This idea is not ap- 
plicable to the recently discovered noncentrosymmetric 
superconductors, CePtaSi, Ulr, CeRhSia, CelrSia, and 
hi2{VAi_^I't^)B,L2A^^ In systems without inversion 
center, parity violation yields the admixture of the spin 
singlet and triplet pairing states, as pointed out by Edel- 
stein nearly two decades agoZ^ An asymmetric potential 
gradient VF caused by the absence of inversion symme- 
try gives rise to the spin-orbit interaction a{k x W) ■ cr 
which splits the Fermi surface into two pieces and aligns 
electron spins on each Fermi surfaces parallel to fc^? x VF, 
with kp the Fermi momentum. A simple example of 
this situation for the Rashba type spin-orbit interaction 
with W parallel to the z-axis is depicted in FIG[T] 
Note that in FIG.l the spin quantization axis is cho- 
sen along fcp X W. On one of the spin-orbit splitted 
Fermi surfaces, say the (-l-)-band in FIG[I] the Cooper 
pair between electrons with momentum k, spin "f and 
momentum — fc, spin | is formed. This state, denoted as 
\k t) I — i), is not a spin singlet state, because the coun- 
terpart of this state |fc i)| — A; t) is formed on another 
Fermi surface and thus the superposition between these 
two states is not possible. Actually, the pairing state 
|fc T)l ~ i) is tbe admixture of spin singlet and triplet 
states as easily verified by, 



|fcT>|-fci> 



{\k])\-ki)-\ki)\-k])) 
{\k])\-ki) + \ki)\-k])). 



The first and second terms of the right-hand side ex- 
press, respectively, a spin singlet state and a spin triplet 
state with the in-plane spin projection S'inpiane equal to 
0. Since we take the spin quantization axis parallel to the 



(-)-band 



(+)-band 




k',a' 
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FIG. f: An example of two-dimensional Fermi surfaces split- 
ted by the Rashba spin-orbit interaction depicted on the k^- 
ky plane. The inversion symmetry is broken along the z- 
direction. The gray arrows represent electron spins on the 
Fermi surfaces for the (-l-)-band and the ( — )-band. The 
Cooper pairs between electrons with momentum k, spin cr and 
momentum —k, spin —a are formed on each Fermi surfaces. 



xy-plane, this triplet state corresponds to the = ±1 
state for the spin quantization axis along the z-direction. 
(This means that the d-vector of the triplet component 
is parallel to the plane.) The above explanation is also 
applicable to general cases with more complicated form 
of W. This unique superconducting state exhibits vari- 
ous interesting electromagnetic properties as extensively 
argued by many authors 

In addition to the superconducting state, the absence 
of inversion symmetry also affects properties of the nor- 
mal state in a drastic way. One remarkable feature ap- 
pears in paramagnetic effects; i.e. there is a van-Vleck- 
like spin susceptibility due to the transition between the 
spin-orbit splitted Fermi surfaces F^iSii^ This effect should 
be distinguished from the usual van Vleck term of the 
orbital susceptibility by the fact that in contrast to the 
usual van Vleck term, the van-Vleck-like susceptibility 
in inversion-symmetry-broken systems varies as a func- 
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tion of temperature like the Pauli susceptibility. The 
existence of the van-Vlcck-like term is crucial for experi- 
mental investigations on superconducting states because 
it exhibits no change even below the superconducting 
transition temperature for both the spin-singlet and 
triplet states when the spin-orbit splitting is sufficiently 
larger than the superconducting gap.^'^ For the deter- 
mination of the symmetry of the Cooper pair, we need 
precise knowledge about the behaviors of the van-Vleck- 
like susceptibility both below and above Tj,. 

Furthermore, the parity-breaking spin-orbit interac- 
tion gives rise to nontrivial coupling between charge 
and spin degrees of freedom, and thus causes unique 
transport properties. For example, there exist mag- 
netoelectric effects such as the electric-field-induced 
magnetizatiouf^i^^ and the spin Hall effect which is char- 
acterized by the transverse spin current caused by an 
electric field.--'?- Also, the anomalous Hall effect consid- 
ered by Karplus and Luttinger five decades ago is possible 
to occuij^^ For heavy fermion systems without inversion 
symmetry, the experimental studies on these effects are 
important future issues. 

In the heavy fermion superconductors CePtaSi, 
CeRhSia, Ulr, and CelrSis, strong electron correlation 
which yields an extremely enhanced effective electron 
mass plays a crucial role. For the clarification of the 
nature of these systems, it is important to understand 
effects of the electron correlation on the aforementioned 
exotic phenomena. The main part of this paper is de- 
voted to this purpose. On the basis of a general frame- 
work of the Fermi liquid theory, we derive the formulae 
for the spin susceptibility and the transport coefficients 
related to the anomalous Hall effect, the thermal anoma- 
lous Hall effect, the spin Hall effect and the magnetoelec- 
tric effect taking into account strong electron correlation 
effects in a formally exact way. Using these formulae, we 
argue how strong electron correlation inherent in heavy 
fermion systems affect experimentally observable quanti- 
ties. The insights obtained from this analysis are useful 
for the exploration of the superconducting state realized 
in CePtaSi, CeRhSia, and CelrSia, as explained below. 

The main results of this paper are summarized as fol- 
lows: 

(i) The thermal anomalous Hall effect, which is the 
anomalous Hall effect for the heat current is predom- 
inated by the contributions from electrons away from 
the Fermi surface. It implies that when the spin-orbit 
splitting is much larger than the superconducting gap, 
the behavior of the thermal anomalous Hall conductiv- 
ity K^j^^ in the superconducting state is similar to that 
in the normal state. Even in the zero temperature limit 
T ^ 0, k/^™" /THz with an applied magnetic field is 
finite, and much enhanced by electron correlation effects 
in heavy fermion systems. 

(ii) In the Rashba case where the inversion symmetry 
is broken along the z-direction, as in the case of CePtaSi, 
CeRhSia, and CelrSia, the spin Hall conductivity, which 
is given by the correlation function between a charge cur- 



rent and a transverse spin current, is not at all affected 
by electron correlation effects, but determined only by 
the band structure. 

(iii) We clarify electron correlation effects on the mag- 
netoelectric effect both in the normal state and in the 
superconducting state. As was considered by Levitov 
et al.^ in the normal state, the bulk magnetization is 
induced by an appfied electric field, M = TE, and, 
conversely, the charge current flow is induced by an 
AC magnetic field, J = —TB. We find that in the 
temperature region where the resistivity exhibits T^- 
dependence, the magnetoelectric effect coefficient T be- 
haves like ^ l/(7r^), with 7 the specific heat coefficient, 
while in the temperature region where the resistivity is 
dominated by impurity scattering, T is enhanced by a 
factor ~ 7. In the superconducting state, the paramag- 
netic supercurrent is generated by the Zeeman field, as 
extensively argued by several authors j^iiiii^iiii We reveal 
electron correlation effects on this supercurrent from a 
general theoretical framework. 

(iv) The temperature dependence of the van- Vleck- like 
spin susceptibility x^^{T) crucially depends on the de- 
tails of the electronic structure, and occasionally may 
become stronger than that of the Pauli susceptibility 
^Pauh^j^-j^ in contrast to the usual van Vleck orbital sus- 
ceptibility, which is independent of temperatures. Also, 
in the superconducting state, the ratio of x^^C^) to the 
total spin susceptibility xiT) = x^''""(T') + X^^ {T) in 
the zero temperature limit is seriously affected by both 
the electronic structure and electron correlation effects. 
This result yields an important implication for the recent 
NMR measurements for CePtaSi carried out by Yogi et 
aZ.,^^'^^ which indicate no change of the Knight shift be- 
low Tc for any directions of an applied magnetic field, 
contrary to the previous theoretical prediction that for a 
magnetic field along the z-axis, x(0) = x(Tc)/2 for a suf- 
ficiently strong spin-orbit coupling. Our precise analysis 
indicates that for a specific electronic structure, x(^) is 
dominated by x^^j ^iid as a result, x(0) ~ xC^c), which 
is consistent with the aforementioned NMR experiment. 

The organization of this paper is as follows. In the 
next section, we develop the analysis for the normal state, 
exploring the static magnetic properties and the trans- 
port coefficients. In the section III, we extend our ap- 
proach to the superconducting with particular emphasis 
on the magnetism. We propose a possible explanation for 
the recent NMR experimental data. We also argue elec- 
tron correlation effects on the paramagnetic supercurrent 
caused by the inversion-symmctry-breaking spin-orbit in- 
teraction. A summary is given in the section IV. 



II. NORMAL STATE 

A. Basic equations 

In this section, we consider a single band interacting 
electron model with a lattice structure breaking inversion 
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symmetry. An extension to more complicated systems 
such as the periodic Anderson model is straightforward, 
and will be discussed elsewhere. Our model Hamiltonian 
is given by, 



H = Hq + Hso, 



k.a 



Ck, 



(1) 
(2) 

(3) 



where — (c|j.,c|j.) is the two-component spinor 
field for an electron with spin |, j, and momentum 
k. n„i = c^aiCai IS thc number density operator at 
the site i. a — {cTx, '^yi '^z) with a^,, v = x^y,z, the 
Pauli matrix. For simplicity, we have assumed the on- 
site Coulomb repulsion with the coupling constant U in 
Eq.(l2]). The vector Co{k) — {Cqx , Coy , Cqz) in the spin- 
orbit interaction term Hso obeys Co{—k) = —Co{k), 
of which the explicit expression is determined by the 
detail of the crystal structure which breaks the inver- 
sion symmetry. For the tetragonal symmetry and small 
k, Co{k) = {ky,—kx,0) which leads the Rashba-type 
spin-orbit interactioni^ For the cubic symmetry with 
Zinc Blende structures and small k, Co{k) = {kj;{ky — 
k1), ky{k'^ — kl), kz{kl — ky)), which is related to the Dres- 
selhaus spin-orbit interactioni^i^i^ 

In the following, we are concerned with magnetic prop- 
erties of the system. Unique features of magnetism in 
inversion-symmetry-broken systems appear as a result 
of paramagnetic effects rather than diamagnetic effects. 
Thus we add to the Hamiltonian ([1]) the Zeeman coupling 
with an external magnetic field H — {Hx, Hy, Hz), 



-H^Zccman — ^ C^.^B"' ' Hck, 



(4) 



and neglect the orbital diamagnetic effects. For simplic- 
ity, we assume that the g-value is equal to 2. Then, the 
inverse of the single-particle Green's function under the 
applied magnetic field is defined as. 



G ^(p) = ie„ - H{p) + /iBcr • H, 



where p = {iSn, k), and, 



i/(p)-i/o(p) + S(p), 



(5) 



(6) 



Here S = (S„I]y,S,) with Sq = (E^ + Su)/2, 
(E^t + STi)/2, E, = (E^T - 5^Ti)/2*, and E, = (Ett 
Sii)/2. 

It is convenient to introduce a vectorial function, 



C{p) = {Cx{p),Cy{p),Cz{v)) = ^o(fc) - + -S(p).(9) 



a 



Generally, C{p) is non-Hermitian, because E^; and E^ 
may be complex quantities. However, G^^{p) = ie — £fe + 
/i— Eo(p)— q:/^(p)-c is diagonalized by the transformation 
Uip)G-\p)U+ip) with. 



Uip) 



u+ip) 



-^-{p)v+ip) C+b) 



^+{p) -^£.-{p)v-{p) 
i-{p)v+{p) ^+{p) 



1 ± 



\\C{p)\\ 



■n±ip) 



^ Cx{pf + Cy{pf 



(10) 



(11) 



(12) 



(13) 



Here ||/^(p)|| = y '^^ + ^y+ ^1- the absence of the 
magnetic field, we can verify that the time reversal sym- 
metry leads the relation 



E* (z,-fc) = -E^t(^>fc)- 



(14) 



As seen from eqs.®,©, and (jS]), the main effect of 
the off-diagonal self-energy terms, E||, E|| is to renor- 
malize the spin-orbit interaction term. Since the on-site 
Coulomb interaction does not change the symmetry of 
the system, the off-diagonal self-energy should obey the 
same symmetric properties as Coik) in the absence of 
magnetic fields; i.e. 



S(ie„,-fc) = -S(i£„,fc). 
In particular, for the Rashba interaction, 

Ex(fca;j ky) — '^^xikxjky), 



(15) 



(16) 



^0 = £fe ~ P + aCo{k) ■ a (7) 

with [i the chemical potential. The self-energy matrix E 
consists of both diagonal and off-diagonal components, 

- _ / Eti(p) E|x(p) \ 
^ - \ Y.,^{p) Eu(p) j 

= Eo + S -cr. (8) 



Ey( kx,ky) — Yiy(^kx,ky), (f^) 

and Ej, (Ej^) is an even function of kx (ky). 

The relations and imply that when thc time 
reversal symmetry is preserved, E^, and E^ are real quan- 
tities, and thus U{p) is unitary, i.e. U+{p) — UHp)- In 
more general cases with H 0, as long as the spin-orbit 
splitting is much smaller than the Fermi energy, as in 
the case of any heavy fermion systems without inversion 
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symmetry, the off-diagonal terms of E may be negligi- 
ble compared to the diagonal terms, and [/+ is safely 
approximated by . 

The single-particle excitation energy e^.^ for the quasi- 
particle with the helicity t = ±1 is given by the solution 
of the equation Det[G'~^(£^^, fc)] — 0, which is, in the 
diagonalized representation, 

el, -f M - e/c - So(4, k) - ™||£(4„ fc)|| = 0. (18) 

Using the transformation U (p) , we can easily obtain the 
inverse of Eq.®, 



G{p) = 



1 + r£(p) • a 



T = ±l 



Grip), 



where C{p) — C{p)/\\C{p)\\, and 

Grip) ^ 



(19) 



(20) 



fi-T,oip) - Ta\\Cip)\\ 

In the vicinity of the Fermi surface, the quasiparticle 
approximation is applicable, and the Green function for 
= is reduced to. 



Gip) = J2 



l + r£(£*.,,fc).^ 



Grip), 



r=±l 



with 



Grip) 



Zkr 



i£n + «7feTSgne„ - £ 



(21) 



(22) 



kr 



Here ^kr is the quasi-particle damping given by 



Ikr 



ZkrilmT,^ie, k) + aReC ■ ImS^(£, fc))|e^,e„(23) 



where KeC = Re/^^(e, fc)/|Re/^^(£, fc)| and quantities 
with the superscript R indicates retarded functions ob- 
tained by the analytic continuation ie„ — s- e + i(5 (5 > 0). 
The mass renormalization factor Zkr is 



Zkr 



1 



aso(p) d\\cip)\\ 

— Ta- 



diiSn) 



d{i£n) 



(24) 



In particular, in the case of the Rashba interaction with 
a potential gradient along n = (0, 0, 1), 



Zkr 



t. 



^ dJ^oip) _A as.(p) 

d{ien) \ d{ien) 
dT.yip) 
" d{ien) 



(25) 



with a two-dimensional vector tip) — t(p)/p(p)||, f(p) = 
n X Cip), and \\tip)\\ = \Jt'i+ ty- In the derivation of 

Eqs. (|23)) . ((24)l . and ((25)) . it is assumed that the quasipar- 
ticle damping is sufficiently small to stabilize the Fermi 
liquid state; i.e. Iml]^(e, A;) ^ e with = 0,x,y,z. We 
would like to stress that the notations Ti^^y is more con- 
venient than the notations and because ImS|| 
and IniS|':|- are not directly related to the quasiparticle 
damping, and can not be assumed to be small. 



B. Ward identities 

In the following sections, we will argue transport phe- 
nomena and magnetic properties. In the derivation of 
transport coefficients, the Ward identities which relate 
three-point vertices of multipoint correlation functions 
to the single-particle Green function play a crucial role. 
These identity relations are derived from the conserva- 
tion laws for charge and spin degrees of freedom. 

The Ward identity corresponding to the charge conser- 
vation law has the standard form, 

J2 q^Alip + q,p)^G-\p)-G-\p + q), (26) 
^—0,x,y,z 

where = i—iuj, q), and each components of are the 
three-point vertex functions dressed by electron-electron 
interaction for the charge density (/i = 0) and the charge 
current (/i — x, y, z) In the limit of g ^ 0, the charge 
current vertex function is given by. 



d 



A^(p,p) = ^i£k + So + oiicT ■ Cip))). 



(27) 



In contrast to Eq. ([26]) . the Ward identity for the spin 
degrees of freedom is affected by the spin-orbit interac- 
tion which breaks the conservation of the local spin den- 
sity. To derive the Ward identity for the spin degrees of 
freedom, we start from the continuity equation for the 
spin density. 



dS^jx) 
dt 



v-r''ix)^Pl{x). 



(28) 



Here, in terms of the spinor field ijj'^ix) ~ 

cj, exp(— ifca;), the z^-component of the spin 

density operator S^ix) [v — x, y, z) and the correspond- 
ing spin current density operator J'^'^ix) are expressed 



S-ix)^i;\x)^ijix), 



(29) 



J"" = -[i>Hx)<y.^i'ix) - {v^t(^)}^^^(^)]. (30) 

4i 

In (j28|) . Pgix) is the source term caused by the spin-orbit 
interaction Tiso- From the coordinate representation of 

/dx 
— [ijHx)cT ■ {Coi^iV)^ix)} 

- {Coi-i^)i^\x)} ■ cT^ix)], (31) 
we have the explicit form of Pg(a;), 

Psi^) = ^[^Hx)(T ■ in, X Coi-tV))^jix) 

~ {in, X Coi-iV)) ■ (7TpHx)}^Pix) 

- il^\x)iCQyi~iV)%ljix) 

- i{Co,i-iV)i^\x)}^ix)l (32) 
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where is a unit vector along the i/-axis. In the case 
of the Rashba spin-orbit interaction which breaks the 
inversion symmetry in the (0, 0, l)-direction, /^o(fc) = 
to(fc) X riz- Here the vector to{k) = {tox,toy,0) trans- 
forms Uke {kx, ky,0). Then, the source term for is of 
the form, 



-{to(-*V)^t(^)}.^^(^)]^ (33) 



and for S*^, 



P2 ''2 



ai02 (Pl>P2;P3.P4) = 



Pi Ol 



Pa 03 



P4 04 



FIG. 2: Diagrammatic representation of the four-point vertex 
with spin indices ai, (T2, 0-3, and 0-4, which is irreducible with 
respect to particle-hole pairs. 



The standard argument applied to the continuity equa- 
tion leads the Ward identity for the spin degrees of 
freedomi^^ 



fi—0,x,y,z 



2 
(35) 



where A^'^(p -I- q,p) {fi = 0,x,y,z) are the three-point 
vertex functions dressed by electron-electron interaction 
for the i/-component of the spin density (/i — 0) and the 
corresponding spin current density (/i = x, y, z). T'^{p + 
q,p) is the dressed three-point vertex function for the 
spin source pg. A^'^ and T'^ are the 2x2 matrices in the 

spin space. The three-point vertex T'^{p + q,p) is given 
by the solution of the following integral equation. 

f'^{p+l.p-\)=%^{p+\,p-\) 
p' 

xG{p' + \)f'^{p' + - |)G(p' - |)]. (36) 

Here ^(pi,P2]P3iPi) is the four-point vertex function 
which is irreducible with respect to particle-hole pairs. 
The diagrammatic representation of T with spin indices 
is shown in FIG [2] The bare vertex function is defined 
as 

%%P + I'P - f ) - f ['Co.lfc - f ) - ^o.(fe + f )] 
+ ^({£o(fc + f ) + C,{k - |)} X a) ■ n,. (37) 

In the case of the Rashba spin-orbit interaction, the z- 
component of the bare vertex function is, 



ro^(fc-f |,fc-|) = ^^.[t„(fe + |) 

and the x-component is. 



to(fc-|)],(38) 



%^{k + |, fc - |) = -|[t0,(fc + f ) - hy{k - |) 



(39) 



The above formulae (|55|) — (|59|) are useful for the compu- 
tations of various correlation functions considered in the 
following sections. 

In the limit of 5 ^ 0, the Ward identity ([55)) is reduced 
to the following relation, 



T''{p,p) = ia{C{p) X cr) • n^. 



(40) 



The above relation is not so trivial. Thus it is worth while 
examining (|40p in more details. Substituting Eq. (|l(I)l into 
Eq. (|36[) with (7 = 0, and using the following identity, 
which is derived from (fT9|l . 



G(p)*a(£(p) X <T),G(p) = _^ 



Gnip) 

-uip) 



(41) 



we find that the identity (j40p with v = z is equivalent to 
the following relations for the self-energy functions, 

Stx(p) - Y.lT\l{p,p;p',p')Gn{p') 



p' 



-T\l{p,p;p',p')G^^{p% (42) 

SiT(p) = Y.K\(P'P--P''P')Gll{p') 

p' 

- rl]{p,p;p',p')Gn{p% (43) 

J2^\]{p,p;p',p')Gn{p') = 5]r||(p,p;p',p')GiT(pl,44) 
p' p' 

5]r|I(p,p;p',y)GTi(p') = Eru(P:P;p'.p')GiT(p'I-45) 



p' 



p' 



Here F^iJ^J is four-point vertex functions F which ap- 
pear in eq. (l36p . The interpretation of the identity re- 
lations in terms of the diagrammatic language 
is as follows. In FIGlSl we depict diagrams which rep- 
resent the off-diagonal self-energy S-f^. It is noted that 
the four-point vertex functions F in FIGl2]are obtained 
by rotating the four-point vertex shown in FIGl2] by 90 
°. Thus, the four-point vertex functions in FIG 13] are 
not irreducible with respect to particle-hole pairs in the 
horizontal direction. Among these diagrams, FIGsj^I^a) 
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(a) 



(0 



(b) 



id) 



(a) 



(b) 



(c) 



1 , „ P+q/2| i p'+q/2 P+q/2 I f p'+q/2 

k.p-pjQk 




(g) 



(h) 



FIG. 3: Diagrams representing the ofF-diagonal self-energy 
The gray box is a four point vertex. The up and down 
arrows, respectively, express the up and down spins. 



and are also expressed in the form of FIGlSlJc), 

since in the four-point vertex in FIGsO^a) and ^h), 
there should be at least one G|j-line. We note that 
the way of expressing in the form of FIGlHIc) is 
not unique. The double counting occurs for diagrams 
which are also rewritten in the form of FIGl2)^d), be- 
cause the four point vertex in FIGlH^d) contains two 
more G^j^ -lines than that in FIGl31[c), and thus the di- 
agrams of the form FIGlSfd) are expressed in the form 
of FIGini^c) in two different ways. We show in FIG|4] 
some examples of this double counting which occurs for 
Ep.r||(p,p;p',p')GTi(p') etc. In FIGsffla), (b), and 
(c), we depict some examples of the second-order dia- 
grams for T'^ {p + q/2, p~ q/2). The Ward identity relates 
these three different diagrams, respectively, to the dia- 
grams FIGsllJd), (e), and (f), which express the second- 
order terms of S^j. The diagrams FIGsHI^d) and (e) 
are parts of the diagram FIGlSjc), and the diagrams 
FIGg);f) is a part of the diagram FIG^d). Appar- 
ently, these three self-energy diagrams FIGslUJd), (c), 
and (f) are equivalent to each other, and double count- 
ing of diagrams occurs in FIGO^c). Thus, S|| is repre- 
sented by FIGlSl^c) from which the duplicate diagrams 
of the form FIGl3]Jd) are subtracted. This amounts to 
the identity We can also give a similar interpre- 

tation to Eq.(I321)- The imphcation of Eqs.dMI) and (liS)) 
is more transparent. The left- and right-hand sides of 
Eq. (|44p are nothing but two equivalent representations of 
In FIGHl wc exhibit some examples of this equiv- 
alence. FIGsHJg) and (h) are examples of diagrams for 
T'^{p + q/2,p — q/2), which are, respectively, related to 
the self-energy diagrams FIGsHJi) and (j), which con- 
stitute These equivalent self-energy diagrams are, 

respectively, a part of J2p' ^\]{PtP'jP' jP')^uip') a 
part of J2p' r||(p,p;p',p')G|-f (p'). This amounts to the 
relation The relation holds for any higher-order 

diagrams. Eq. ([^51) also expresses I^n in two different 
manners. 



p+q/2 f ^ i p'+q/2 P+ q/2 f I p'+q/2 

I \^ I p -q/2 



P+P'-k 

p-q/2 i f 

(i) 




p'-q/2 



SO' 

•or 



k+p-p' 

p-qo 
(j) 



FIG. 4: (a), (b), (c) Examples of second order diagrams 
which constitute T'^{p + q/2,p — q/2). Broken lines repre- 
sent the Coulomb interaction U. Black triangle represents 
the bare vertex Tq . (d), (e) Second-order self-energy dia- 
grams which constitute ^ , r|J(p,p;p',p')G|x(p'). The dia- 
grams (d) and (e) are, respectively, related to the diagrams 
(a) and (b) via the Ward identity, (f) A second-order self- 
energy diagrams which constitute ^ , r^^^{p, p; p' , p')G n {p') . 
The diagram (f) is related to the diagram (c) via the Ward 
identity. Note that the self-energy diagrams (d), (e), and 
(f) are equivalent to each other, (g), (h) Examples of third- 
order diagrams for T"^ {p + q/2,p — q/2). (i) A third-order self- 
energy diagram which constitute rj|(p,p;p', j3')G||(p'). 
The diagram (i) is related to the diagram (g) via the 
Ward identity, (j) A third-order diagram which constitute 
, r||(p,p;p',p')Gtx(p'). The diagram (j) is related to the 
diagram (h) via the Ward identity. Note that the diagram (i) 
is equivalent to the diagram (j). 



C. Magnetism: Pauli and van-Vleck-like spin 
susceptibilities 

As mentioned in the introduction, one of the unique 
magnetic properties in inversion-symmetry-broken sys- 
tems is the existence of a van-Vleck-like spin susceptibil- 
ity which stems from magnetic response of electrons oc- 
cupying the momentum space sandwiched between spin- 
orbit splitted two Fermi surfaces. In this section we ob- 
tain general formulae for the Pauli and van-Vleck-like 
spin susceptibilities taking into account electron correla- 
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tion effects in a formally exact way. The uniform spin 
susceptibility is easily obtained from p9p and (PO]) by 
using 



as,, 



(55) 



-rEE^Btr[a^G(p)]|H, 



(46) 



e„ A: 



In the case of cubic systems, a straightforward calcu- 
lation gives, 



Mb 



EE 



-4Tcosh2(^) 



— -— r— Ap"'^(e^^, fe) 



Mb 



EE-^|Re£(e*^,fc)| 



A^"'^(4,,fc),(47) 



with /(e) the Fermi distribution function. Here the 
three-point vertex functions A™^ and Ay"*^ are. 



1 9Eo(p). 



MB ai/^ 

/IB dH:, fiB dH^ 



(48) 



{Cl+Cl){\- 



1 9S]o(p) 



/^B ai72 



MB 



Mb 



OH, 



(49) 



The first and second terms of (|47|) are, respectively, the 
Pauli and van-Vleck-like terms. The van-Vleck-like term 
is not much affected by superconducting transition if 
the magnitude of the spin-orbit splitting Q;|/^(p)| is suffi- 
ciently larger than the superconducting gap. 

In the case of the tetragonal symmetry with a poten- 
tial gradient along n = (0,0,1) (i.e. the Rashba inter- 
action), Co{k) is parametrized as to{k) x n. The vector 
itox,toy,0) transforms like {kx, ky,0). Then, using a vec- 
tor t{p) — n X C{p) and a unit vector t{p) = t(p) /\t{p)\, 
we express the spin susceptibility as. 



Xz 



T = ± k 



a|Re *(4„fc)| 



A-(4„fc), (50) 



Pauli 

Xxx 



_ Pauli , VV 
Xxx — Axx ^ Xxx 1 



' = ^bEE — (52) 



(51) 



■,^V4rcosh2(^) 



,VV _ ,,2 ^ ^ . '^'^^^^kr) 



Xxx — Mb 



(4,,fc). (53) 



The three-point vertices A'^^ , A^^_ , and A^^ in the above 
expressions are given by 



A-(p) = 1 



1 d^ojp) 
MB dHz 



(54) 



A-_(|,)-L(l 



1 d^x. 
p.B dHx ' p.B dHx 



ty 9S.y 



(56) 



As seen from (|50|) and ((5T|) , Xzz is given only by the van- 
Vleck-like susceptibility, while Xxx consists of both the 
Pauli and van-Vleck-like terms. 

We note that in the limit of a — > (the inversion sym- 
metry and the spin rotation symmetry are recovered). 



d^x 



5So 
OH' 



(57) 



This relation is easily verified by rotating the spin axes 
as (Ta; — > fJz, — > — fa;, and (Jj, — > (Jj,, and noticing that 
this rotation transforms (S||-|-S||)/2 into (S||— E||)/2. 
Thus, Eqs. ((50|) and (jSTj) leads Xzz = Xxx in the limit of 

We see from ([50)1 and ([^T]) that both the Pauli term and 
the van-Vleck-like term are enhanced by the factors A*^ , 
A^^_, and A^^. If the density of states weakly depends 
on energy, and the spin-orbit splitting is much smaller 
than the Fermi energy, these three enhancement factors 
take values of the same order. However, when the en- 
ergy dependence of the density of states is substantial, 
the enhancement of the van-Vleck-like susceptibility due 
to the electron correlation effect may be quite distinct 
from that of the Pauli term. We will consider such an 
example in Sec. HI in connection with the heavy fermion 
superconductor CePtsSi. 

For experimental studies on superconducting states, it 
is important to discriminate between the Pauli suscepti- 
bility which decreases in the spin-singlet superconducting 
state, and the orbital susceptibility which is not affected 
by the superconducting transition. Usually, this discrim- 
ination is achieved by experimentally observing the tem- 
perature dependence of the spin susceptibility in the nor- 
mal state, and identifying the temperature-independent 
part with the orbital term. However, this idea fails for 
noncentrosymmetric superconductors because of the exis- 
tence of the van-Vleck-like susceptibility x^J- When the 
spin-orbit splitting is much larger than the superconduct- 
ing gap as in the case of CcPtaSi, CcRhSia, and CelrSis, 
the van-Vleck-like term x^x ^'^^ much influenced by 
the superconducting transition and takes a finite value 
even at T = 0, though it exhibits temperature depen- 
dence at high temperatures like the Pauli susceptibility. 

Here we demonstrate to what extent x^x varies as a 
function of temperature by using a simple model. We 
consider a two-dimensional electron system on a square 
lattice with the energy momentum dispersion relation. 



Ek — — 2(cOsfcx -I- COS ky), 



(58) 



and the Rashba spin-orbit interaction a{tQ{k)xn)-(T. We 
put the hopping integral t = 1 in (j58p . To preserve the 
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periodicity of the energy band in the momentum space, 
we assume the form of to(fc) s-s 

to(fc) = (sin^,sin^,0). (59) 

This model can not be directly related with heavy 
fermion superconductors such as CePtaSi, CeRhSia, etc. 
which have the complicated three-dimensional Fermi sur- 
faces. Nevertheless, this simple model is quite instructive 
in that it demonstrates the strong temperature depen- 
dence of J ■ For simplicity we neglect electron-electron 
interaction, and consider the half- filling case n = 1, in 
which the van Hove singularity of the density of states 
gives rise to the temperature-dependence of the spin sus- 
ceptibility. For a 7^ 0, to keep the electron density n — 1, 
we need to add a counter term 

-^/i^4cfe, (60) 

k 

with Sfi = — to the Hamiltonian. The calculated re- 
sults of x^™'' and x^J for some values of a are plotted as 
a function of temperature in FIGlSl It is seen that as T 
increases, the curves of xjfj and x^™'' approach to each 
other, and for T 3> a, there is no substantial difference 
between them. The van-Vleck-like term x^:^ varies as 
a function of temperature, reflecting the singular energy 
dependence of the density of states in the model ([55)1 . 
If we take into account electron correlation effects, the 
energy dependence of the density of states is much en- 
hanced, and thus, as T is changed, xS' varies more dras- 
tically than the non-interacting case. This calculation 
suggests that it is almost impossible to distinguish the 
van-Vleck-like term from the Pauli susceptibility merely 
by analyzing experimental data of temperature depen- 
dence of the spin susceptibility in the normal state. 

It is also worth while examining the anisotropy of the 
spin susceptibility due to the spin-orbit interaction. The 
calculated results of Xxx and Xzz as a function of T are 
plotted in FIGlH] for several values of a. It shows that 
the anisotropy becomes large for T < a|tofe|- The mag- 
nitude of the anisotropy \xxx — Xzz\IXxx is at most of 
order ^ a\tQk\/EF: as easily expected from (|50)) and 
([5T|) . According to the LDA calculations for the heavy 
fermion superconductor CePtsSi, the ratio a\tQk\/Ep is 
of order ~ q ]^, 34^35 Although the de Haas van Alphen 
experiment for this system has not yet detected the spin- 
orbit splitted Fermi surfaces, the measurement for a re- 
lated system LaPtsSi support the existence of the spin- 
orbit splitting of which the magnitude is of this order 
Thus, it is expected that the anisotropy of the spin sus- 
ceptibility for CePtaSi is at most ~ 0.1, which seems to 
be consistent with the experimental observation of the 
susceptibility,^!^ 

D. Transport coefficients 

Generally, the absence of inversion symmetry brings 
about remarkable transport properties such as the 
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FIG. 5: Xiix"'' and Xxx versus T. Each lines are, respectively, 
for a = 0.1 (top solid line), x^^S"'' for a = 0.1 (top 
dotted line), x^^ for ot — 0.5 (bottom solid line), and XiS"'' 
for Q = 0.5 (bottom dotted line). 
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FIG. 6: Xxx and Xzz versus T. Each lines are, respectively, 
Xzz for a = 0.1 (top dotted line), Xxx for q = 0.1 (top solid 
line), Xzz for a = 0.5 (middle dotted line), Xxx for a = 0.5 
(middle solid line), Xzz for a = 1.0 (bottom dotted line), and 
Xxx for a = 1.0 (bottom solid line). 

anomalous Hall effect ^i^^ the thermal anomalous Hall 
effect, the spin Hall effect r^i^l and magnetoelectric 
effects. ^Ai^ In this section, we argue these phenomena 
with particular emphasis on the role of electron correla- 
tion effects. We mainly consider the case of the Rashba 
spin-orbit interaction with Co = x n in which the 
above-mentioned effects are more salient than the case of 
cubic systems with the Dresselhaus interaction. 



1. Anomalous Hall effect 

Here, we consider the anomalous Hall effect due to the 
Rashba spin-orbit interaction. The general theoretical 
description of the anomalous Hall effect due to the spin- 
orbit interaction was given by Karplus-Luttinger many 
years ago.^^ The key ingredient of the Karplus-Luttinger- 
type anomalous Hall effect is the existence of a dissi- 
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pationless transverse current which is induced by the 
paramagnetic Zeeman effect and stems from the anoma- 
lous velocity, Vfea(to(fc) x n) • cr. In the case of the 
Rashba spin-orbit interaction, there is no z-component 
of the anomalous velocity, and thus this effect exists 
only for an applied magnetic field along the z-dircction. 
In heavy fermion systems, the spin-orbit interaction 
caused by heavy atoms like Ce and U may also give 
rise to the anomalous Hall effect as experimentally ob- 
served in several systems such as CeRu2Si2, CeAls, and 
UPta i^'^i^'^i'^ The strong anisotropic dependence on 
an external magnetic field of the Rashba-interaction- 
induced anomalous Hall effect makes a clear difference 
from the conventional anomalous Hall effect observed in 
the systems with inversion symmetry. Note that in the 
case of cubic systems with the Dresselhaus interaction, 
such a strong anisotropy does not exist. 

The Hall conductivity for a magnetic field in the z- 
direction is defined as 



7xy — lim —K.j^y(iu}n)\iu]„^uj-\-i<ii 



(61) 



corrections are obtained by extracting the particle-hole 
pan G^G^ with G^(^) the retarded (advanced) single- 
particle Green function after the analytic continuation 
of Kxy{iio,i), since only G^G^ is singular in the limit of 
bj for infinitesimally small quasiparticle damping. — 
In our case, as will be shown below, the particle-hole 
pair G+G+ and G-G- do not enter into the expres- 
sion of the anomalous Hall conductivity, and only the 
pair G+G^ contributes. When the spin-orbit-splitting 
is sufficiently larger than the quasiparticle damping, the 
particle-hole pairs G^G^ and G^G^ are not singular, 
and do not yield dissipative processes. Thus, we can ig- 
nore the dissipative vertex corrections for systems with 
large spin-orbit splitting. As a matter of fact, for the non- 
centrosymmetric heavy fermion superconductors CePtsSi 
and CeRhSis, the typical size of the spin-orbit-splitting 
is of order 0.1 Ep with Ep the Fermi energy, and dom- 
inates over the quasiparticle dampingi^i^i^ Then, the 
current-current correlation function K^y is expressed in 
terms of the Green functions. 



Kxy{iuJn) 



l/T 



dT{Tr{MT)Jy{0)})e'^-\ (62) 



The total charge current is, 



where the velocity Vkf_i is 



Vkf, = Vfe^(efc + acr ■ {to{k) x n)). 



(63) 



(64) 



Since the anomalous Hall effect is caused by the in- 
terplay between the spin-orbit interaction and the Zee- 
man coupling with an external magnetic field, we ne- 
glect the orbital diamagnetic effect of the magnetic field 
which yields the normal Hall effect. Then the Hamilto- 
nian under consideration is given by Eqs.HI) and (g]) with 
H — (0, 0, Hz). From Eq. (|?7)l . the charge velocity vertex 
fully dressed by electron-electron interaction is given by. 



A^(p,p) = Vfe^(£fe + So(p)) + Vfc^E,(p)(7^ 



(65) 



The last term of the right-hand side of ((65)) is the anoma- 
lous velocity. In the above expression of Aj^, the vertex 
corrections associated with dissipative processes which 
are important for the total momentum conservation^^ are 
not included. The dissipative vertex corrections cannot 
be taken into account by the procedure presented above. 
However, their contributions to the anomalous Hall con- 
ductivity are negligible in the case that the spin-orbit- 
splitting is much larger than the quasiparticle damp- 
ing because of the following reason. According to the 
general argument of the Fermi liquid theory for non- 
equilibrium transport coefficients, the dissipative vertex 



XVkyG{p)]. 



(66) 



Here q = (io;, 0,0,0). Collecting terms linear in for 
Hz — > 0, we find that only the anomalous velocity term 
in Eqs. (|64p and (j65[) gives non- vanishing contributions to 
Kxy{ioj), which have the form J^p /(fc)tr[cr'^G(p)cr^G'(p-|- 

q)] or Ep/(fc)tr[CT^G(p)cr^G(p-hg)], where /(fe) is the 
even function of momentum k. Using the assumption 
that the quasiparticle damping is much smaller than the 
spin-orbit-splitting, we obtain. 



dH, 



■ fik)ir[ayG{p)o-G{p + g)]|H,^o 



p ^ 

X (G+ {p)G- {p + q)- G- {p)G+ ip + q)) 

we(p +q)- eip)e+ip + q))\H.^o 

x{G+{p)G+{p + q)-G^{p)G^{p + q))]. 



(67) 



The last term of the right-hand side of ([67|) does not 
contribute to Re axy- In the limit of w — s- 0, 



■ (a {p)e+ {p+q)-e- {p)e {p+q))\H. 



dH, 



a|Re*(p)P'^° dHz 
Using (|66|) . (|67)) . and ([68]) . we end up with. 



(68) 



Re a 



AHE 

xy 



Hz 



= e fiB 



EE 



-r/(£*.JA-(e*,,fe) 



^ 2a|Rei(4,,fc)|3 
>^{dk^txTdkytoy — dk^tyrdkytox)- (69) 
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Here dk^Ur = dkjv{e, fc)|£=e*^ ; i-e- dk^ does not operate 
on in the argument of (p) . Since we have postulated 
that the spin-orbit sphtting is much larger than the quasi- 
particle damping, the anomalous Hall conductivity (^^^^ 
is not involved with any relaxation time, and thus deter- 
mined only by dissipationless processes. It is noted that 
the anomalous Hall conductivity is enhanced by the fac- 
tor A'*^ which is equivalent to the enhancement factor of 
Xzz, Eg. ([50]). In the heavy fermion system CePtaSi, this 
enhancement factor is of order ~ 60, and the detection of 
the anomalous Hall effect is feasible in such strongly cor- 
related electron systems. We would like to stress that in 
the expression of the anomalous Hall conductivity 
not only electrons in the vicinity of the Fermi surface but 
also all electrons in the region of the Brillouin zone sand- 
wiched between the spin-orbit-splitted two Fermi surfaces 
contribute, in accordance with the fact that Xzz is dom- 
inated by the van-Vleck-like susceptibility. 



2. Thermal anomalous Hall effect 

We now consider the thermal anomalous Hall effect, 
which is the anomalous Hall effect for the heat current. 
To simplify the following analysis, we assume that the 
energy current due to the interaction between quasipar- 
ticles is negligible, and thus the heat current is mainly 
carried by nearly independent quasiparticles. We would 
like to discuss the validity of this assumption in the end of 
this section. Then, we can obtain the thermal anomalous 
Hall conductivity by using a method similar to that used 
in the previous section. Using the heat current related 
to the single-particle energy, 

^Qm = A\[^k^H{p) + H{p)vk^.]ck, (70) 

k 

we define the Hall conductivity for the heat current as, 

«.. = ^(4^J-E4i^^<°^^4i^)' (71) 

where L^^i, is equal to the conductivity tensor a^y^ and, 
ij:::) = hm -4™)(zc.„)kc.^c.+«, m = 1, 2, (72) 

uj—fO lU) 

4i)(zc.„)- C''^ dT{Tr{JQ^.{T)Mme'^-\ (73) 
Jo 

f'^ dT{TAJQ,{T)JQM}y''"^ ■ (74) 
JO 

Extending the argument in the previous subsection 
straightforwardly to the present case, we find that the 



anomalous Hall effect gives rise to the following terms to 
r(i) „„H r^^) 

Ijxy (AnLL J-Jxy 1 
J (m)AHE 

^ T = ± k 

A''^(e* k) 

''2a|Ret(r*^,fc)|3(^'=^^--^^«*°^ " V<9.„MT5) 

with m = 1,2. Then, the expression for the thermal 
anomalous Hall conductivity k^^^ is given by Eqs. ([7T|) 
and ([751). As in the case of the anomalous Hall effect 
for the charge current, the thermal anomalous Hall con- 
ductivity is also dominated by the contributions from 
electrons occupying the momentum space sandwiched be- 
tween the spin-orbit-splitted Fermi surfaces. This prop- 
erty brings about a remarkable effect in superconducting 
states. In the superconducting state, when vortices are 
pinned in the mixed state, the Hall effect for the charge 
current does not exist. Instead, the thermal Hall effect for 
the heat current carried by the Bogoliubov quasiparticles 
is possible. Since we consider the magnetic field perpen- 
dicular to the xy-plane, electrons in the normal core do 
not contribute to the thermal transport in the direction 
parallel to the plane. Below the superconducting transi- 
tion temperature axy is infinite while L*-^-* is finite. Thus 
the thermal Hall effect is governed by the first term of the 
right-hand side of ([7T|) . i.e. the coefficient Lxy ■ In con- 
trast to the normal Hall effect for the heat current which 
decreases rapidly in the superconducting state, the co- 
efficient Lxy is not affected by the superconducting 
transition when the magnitude of the spin-orbit-splitting 
is much larger than the superconducting gap as in the 
case of CePtsSi and CeRhSia. Thus, even in the limit of 
T ^ 0, K^^^/{H^T) takes a finite value. Moreover in 
heavy fermion systems, the magnitude of Kxy^/i^zT) 
in the limit of T — + is expected to be much enhanced 
by the factor A*^. 

Finally, we discuss the validity of the disregard for the 
heat current carried by the interaction between quasi- 
particles. In the vicinity of the Fermi surface, the quasi- 
particle approximation is applicable, and the interaction 
between quasiparticles is much reduced by the wave func- 
tion renormalization factor zf.^ , and may be negligible for 
heavy fermion systems. However, as seen from Ea.(|75p. 
the thermal anomalous Hall conductivity is dominated 
by the contributions away from the Fermi surface, which 
throws a doubt on this approximation. As a matter of 
fact, the above treatment of the heat current is justified 
as far as the spin-orbit splitting Q;|t(0, fc)| is much smaller 
than the Fermi energy Ep- For Q!|t(0, fe)| ^ Ep, the 
quasiparticle approximation is still valid for all electrons 
in the region sandwiched between the spin-orbit splitted 
two Fermi surfaces. The relation Q;|t(0, fc)| <C Ep holds 
when the carrier density is not so low. This condition is 
satisfied for any noncentrosymmetric heavy fermion su- 
perconductors CePtsSi, CeRhSis, CelrSia, and Ulr. 
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3. Spin Hall effect 

Recently, the existence of the spin HaU effect in the 
Rashba model has been extensively investigated by sev- 
eral authorsi^i^ However, electron correlation effects on 
the spin Hall coefficient has not yet been fully elucidated. 
In this section, we explore the Fermi liquid theory for 
the spin Hall effect. The spin Hall effect is character- 
ized by the transverse spin current induced by an electric 
field. For the Rashba model with the inversion symmetry 
breaking along the z-axis, the in-plane spin current with 
the magnetization in the z-direction is considered. Then, 
the spin Hall conductivity is defined as, 

a^™ = lim — i^^™(^c^„)|^..„^..+^o, (76) 

" 1x1— ►0 

i^SHE(,^„)^ /'^^dr(T.{jr(r)J,(0)})e-"^ (77) 
Jo 

Here the total spin current JJ^ is, 

J^' = ^ E ^U^k.a' + ^'ik.)ck. (78) 

Note that we put the g factor equal to 2. The spin cur- 
rent vertex function fully-dressed by electron-electron in- 
teraction is readily obtained from the Ward identity ([55]) 
which is, in the case of the Rashba interaction, rewritten 
as, 

f-L—0,x,y,z 

{-iuj + Ek+i +So(p+ |) -efe-| - So(p- |))y 

-f • [tip + |) + t(p - |)] + r-(p + |,p - |l.79) 

We note that for ui/q and small q, 'T'^{p + q/2,p — 
q/2) is expanded as, 

'r''{p + l,p~l)^f^p,p) + o{q'). (80) 

The absence of the g-linear term in (|80p is verified as 
follows. For w = 0, in the Feynman diagrams which 
constitute T^(p -I- q/2,p — q/2), the total number of G- 
lines which depend on q/2 is equal to that of G-lines 
which depend on —q/2. Each diagram which is derived 
by differentiating a G{p' + g/2)-line with respect to q 
gives the same contribution as that obtained from the 
differentiation of G{p" — q/2)-\me with respect to —q in 
the limit of g 0. This situation is diagrammatically 
depicted in FlGlTl Thus, all terms of order 0{q) in q ■ 
S/qT^ip + q/2,p — q/2) cancel with each other, which 
leads Eg. ([50]) . Then, putting uj = first and then taking 
the limit of q — > in ((79|) . we have 

A^^(p,p) = MB7^(efc + So(p))a^ (81) 




FIG. 7: Diagrams representing lim^^o VqT''(p-|-g/2,p — g/2). 
The gray quadrangle is the four-point vertex irreducible with 
respect to particle-hole pairs. The black triangle represents 
T" . The gray triangle with a cross is a three-point vertex 
generated by g-derivative. The gray hexagon is a six-point 
vertex irreducible with respect to particle-hole pairs. 

We note that in the expression of A^^, the vertex cor- 
rections associated with dissipative processes are not in- 
cluded, as in the case of the anomalous Hall effect consid- 
ered in Sec.H.D.l. As long as the spin-orbit splitting is 
much larger than the quasiparticle damping, such vertex 
corrections are not important. 

By using ([Ml) and K^'^^{iu!n) is recast into, 

tr[Ar (P + q,p)G{p)vkyG{p + q)]. (82) 

p 

A straightforward calculation leads, 

SHE —Tf{el^)tx{el^,k)vxTdkytoy , ^ 

2«|Ro*(4.,fc)P 

where 

V=5fe.(efe + S?(e,fe))|s=et,^. (84) 

In the expression of the spin Hall conductivity ([55]) . the 
factors VxT and t{p) are not substantially affected by elec- 
tron correlation effects, because the main effect of the 
electron correlation on these factors appears only through 
the deformation of the Fermi surface. Thus, the spin Hall 
conductivity cr^™ is not renormalizcd by electron corre- 
lation effects, but determined solely by the band struc- 
ture. 



4- Magnetoelectnc effect 

The existence of the spin-orbit interaction of the form 
(tofc X \7V) ■ (T gives rise to nontrivial coupling between 
charge and spin degrees of freedom, which results in inter- 
esting magnetoelectric effects. For example, there exists 
the bulk magnetization induced by the charge current 
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flow, as pointed out by Levitov et al. many years agoM- 
To argue this phenomenon, we consider the Rashba-type 
system with VT^ along the z-axis. The magnetization 
along the x-direction is generated by an electric field ap- 
plied in the z/-direction. 



(85) 



The inverse effect is also possible; i.e. an AC-magnetic 
field can yield the charge current flow, 



Jx — 2T2 



dBy 

' dt 



(86) 



The prefactor 2 in ([86)1 is due to the fact that in addition 
to the current directly induced by B, there is a magne- 
tization current c(V x M). It is noted that analogous 
effects in the superconducting state, i.e. the supercur- 
rent flow induced by the Zeeman magnetic field, and its 
inverse effect, are also argued by Edelstein, Yip, and the 
present author.— liiii^iiii 

The formula of the magnetoelectric coefficient Txy in 
the case without electron-electron interaction was ob- 
tained by Levitov et al~>^ Here we consider effects of 
electron correlation on T^^y, which may be important in 
the application to heavy fermion superconductors. In 
the normal state, the magnetoelectric coefficient T^j, 
is a non-equilibrium transport coefficient, to which the 
Fermi liquid corrections are obtained by the method de- 
veloped by Eliashberg for the calculation of the electric 
conductivity.— 

According to the Kubo formula, the magnetoelectric 
coefficient is given by. 



1 



1/T 



dT{TAS^{T)Jy{Q)})e^' 



(87) 



(88) 



Applying the Eliashberg's method,'*'' we can express 
this correlation function in terms of the Green functions 
and the three-point vertex functions for the spin density 
and the charge current density. The spin density vertex 
function required for the calculation of l|88p is obtained 
from the Ward identity (|55|) . Taking the limit of q ^ of 
the Ward identity ([55]) . we have the equation satisfied by 
the fully-dressed three-point vertex function of the spin 
density. 



zU-(p+|,p-|) = ^[zc.-S]o(p+|) 

-f[>c.(p + |)-/:.(p-|)] 

+ ^|({£(p+|) + £(p-|)}x,T).n. 



So(p-|)] 



q 



q 



(89) 



In the above expression, q = (iw,0). Using Ea. (|40|) . we 
extract w-linear terms of the right-hand side of (|89p and 
take the limit w ^ arriving at. 



aEo(p)\ adC,{p) 



2 9(ie„) ) 2 d{ien) 

df-ip + q/2,p~q/2) 



d{iuj) 



(90) 



iu — ^0 



Then, using Eqs. ([5(I)l and (|65p for the spin density and 
charge velocity vertex functions, respectively, and ex- 
tracting the factor G^{p)G^{p) = 2mzU{e - el^)/{uj + 
^ijkr) after the analytical continuation to the real fre- 
quency, we have for small lo, 



4Tcosh2(|^) w + 2tjk. 



with, 



de 2doj 



de 



2duj 



, (91) 



Jyr 



(93) 



lyr = Vyr + T{dkytyriy + dkytxrix) + AX. 



(94) 



Here AJ^ and AX are the contributions from the vertex 
corrections to the charge current which are required for 
preserving the momentum conservation law, and related 
to the quasiparticle damping j^t by the Ward identity 
for the charge degrees of freedom4^ In the following we 
are not concerned with the explicit form of these vertex 
corrections, but merely assume that umklapp processes 
or impurity scattering break the momentum conserva- 
tion, leading to a finite value of the dissipative transport 
coefficient. 

In the expression of Eq. (|9ip . the mass rcnormaliza- 
tion factor Zkr cancels after carrying out the momentum 
sum in the DC limit — > 0, and thus electron cor- 
relation effects appear only through H^t- and the factor 
Ikr/zkT = ImEo(fe, EfeT-). This factor is related to the re- 
sistivity p ^ ImEo ^ cq + AT'^, where the constant cq is a 
residual resistivity due to impurity scattering. In heavy 
fermion systems, the coefficient of the T-square term of 
the resistivity A is almost proportional to the square of 
the specific heat coefficient 7^ , which are much enhanced 
by electron correlation. The factor likr ([92]) is enhanced 
like the specific heat coefficient 7 because of the term 
proportional to —dT,^{p)/de. The other terms in H^t- 
are smaller than this term by a factor a\t{p)\/ Ep <C 1. 
Thus, 



Co 



7 

const. 72T^ 



(95) 
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at low temperatures. In the temperature region where 
the resistivity exhibits T^-dependence, T^y behaves hke 
~ l/(7r^), while in the zero temperature limit T^y is 
enhanced by the factor 7. 



The inverse of Eg. (1961) readily obtained as, 



( G{p) F{p) 



(100) 



III. SUPERCONDUCTING STATE 
A. Basic equations 

The Fermi liquid theory in the normal state developed 
in the previous section can be straightforwardly extended 
to the superconducting state as in the case with inversion 
symmetry.— i^i^ This task for the case of the Rashba 
interaction was done by the author in ref.^^. Here we 
present the formalism for general forms of the spin-orbit 
interaction a/^o(fc) • cr. We add a pairing interaction 
to the model Hamiltonian ([1]), which may be originated 
from higher order interaction processes caused by the on- 
site Coulomb term t/, or any other interactions. 

In the conventional Nambu representation,"^^ the in- 
verse of the single-particle Green's function is defined as. 



where 



H{p) + /iBCT 

-At(p) 



H ^ ^A{p) 

ie„ + H\-p) - Mbc^* 



Here H{p) is given by the same expression as Eq.®. 
However, the self-energy in this case includes both 
the normal Green function Gij)) and the anomalous 
Green function F{p) presented below as internal lines. 

iSn — H{p) + fiBcr ■ H and ie„ -I- H*{—p) — /iec* • 
H in Q^^{p) are diagonalized by the transformation 

Aip)g-'{p)A+ip) with. 



A{p) 



A+ip) 



U{p) 
UXi-p) 

u+{p) 
uH~p) 



(97) 



(98) 



Here U{p) and U^{p) have the same forms as in the nor- 
mal state given by pI7)) and pT|) . respectively. 

To diagonalize the superconducting gap A(p), we need 
the assumption that ImE^, and ImSj, (and thus Im^C^^.y) 
are negligible, which implies U+{p) — U''{p), and is justi- 
fied in the absence of magnetic fields, i.e. H — 0. Under 
this condition, A{p) and A^p) in g~^{p) are diagonal- 
ized by the transformation Ag~^{p)A+ if and only if the 
gap function has the form. 



A{p) = A,{k)i(7y + At{k)C{p) 



(99) 



Eq. ([M)l implies that there is no Cooper pair between elec- 
trons on different Fermi surfaces ek+ and Efc-, of which 
the existence gives rise to pair-breaking effects. 



G{P) = E '-^^^^^G.ipl 



T = ±l 



(101) 



(102) 



T = ±l 



and, 



Grip) 



Fr{p) 



^i-\p) 



Ei+\p)Ei~\p) - Al' 



Ei+\p)Ei~\p) 



(103) 



(104) 



E^^\p) =ieT[ek- ^^ + M±P)+ra\C{±p)\], (105) 

with Akr = As(fc) -I- T\C{k)\At{k). When the quasipar- 
ticle approximation is applicable. 



Grip) = 



Zkrit£ + £lr) 



{ie + ijkrSgne)^ - El^ ' 



Frip) 



ZkrAkr 



{ie + ijkrSgne)^ 



El 



(106) 



(107) 



with Ekr = ^/£*kr + Afcr Akr = ZkrAkr- The mass 
renormalization factor Zkr is also given by the same ex- 
pression as Eq. (|M|) with the self-energy S(p) replaced 
with that in the superconducting state. 

Occasionally, some physical quantities for inversion- 
symmetry-broken systems are governed by electrons oc- 
cupying the momentum space region between the spin- 
orbit splitted Fermi surfaces, as argued in the previous 
sections. In such cases, when the spin-orbit splitting is 
sufficiently large, the quasiparticle approximation cannot 
be generally justified, and it may be required to use the 
complicated expressions Eas. (|103p and p04p . However, 
fortunately, for the noncentrosymmetric heavy fermion 
superconductors CePtaSi and CeRhSia, the spin-orbit 
splitting is much larger than the superconducting gap, 
and thus any physical quantities governed by electrons 
occupying the momentum space region between the spin- 
orbit splitted Fermi surfaces are safely approximated by 
the quantities in the normal state, of which the calcula- 
tion is much easier. 
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B. Paramagnetic properties 

We now argue the paramagnetism of this system on 
the basis of the formahsm presented above. The uniform 
spin susceptibihty is calculated from the i^-component of 
the magnetization 



T 



tr 











5ip) 



(108) 



In the derivation of the spin susceptibility, we ignore the 
magnetic-field-dependence of the superconducting gap, 
because the contributions from dAkr/dHi, to the spin 
susceptibility is vanishingly small when the energy de- 
pendence of the density of states is nearly particle-hole 
symmetric. 

In the case with cubic symmetry, the uniform spin sus- 
ceptibility in the superconducting state is calculated as. 



Xz 



F+{p)F^{p)] 



xA™^(j£„,fe). 



(109) 



The first and second terms of the right-hand side of (I109P 
are, respectively, the Pauli and van- Vleck- like contribu- 
tions. The three-point vertex functions Ap''' and Ay''' 
in Eq. (|109p are slightly different from those appeared 
in the normal state (j48|) and (j49|) because the anoma- 
lous Green function F{p) also contributes to the self- 
energy S(p). However, as long as there is no strong 
spin fluctuation, the corrections to Ap"y due to the ex- 
istence of the superconducting gap are small. In the 
case of |AfcT-| ^ a|Z^(p)|, the quasiparticle approxima- 
tion (|106p and (|107p is not applicable to the van-Vleck- 
like term. Thus, we will not write down its explicit 
form in Ea. (|109p . As a matter of fact, in this case, 
the van Vleck term is well approximated by the formula 
for the normal state (I47p . In the case of a spherical 
Fermi surface without electron-electron interaction, and 
£„(fc) = (fc,(fc2 _ kl),kyikl - kl),k.{kl - kl), if the 
condition |Afe^| <^ a|£(p)| ^ Ep is satisfied, the spin 
susceptibility at zero temperature is Xzz(O) ~ fxzzCTc), 
since the van-Vleck-like term remains finite. 

In the Rashba case, for a magnetic field perpendicular 
to the plane, the uniform spin susceptibility is given by, 

Xzz(T) = -2^ilTY,Y.^G+{p)G-{p) + F+{p)F^{p)] 



xA''^(i£„,fc). 
For an in-plane magnetic field^ 



(110) 



ZkT 



, 4T^osh^^^^^^(--'^) 



2M|r^5][G+(p)G_(p) +^^+(p)i^_(p)] 



x4Ai'L(ie„, fc). 



(Ill) 




FIG. 8; A schematic figure of an example of the density of 
states which is significantly enhanced in the energy region 
sandwiched by the spin-orbit splitted two bands ek+ and Ek-- 
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FIG. 9: Xx'^{T = G)/x^x{T = T^) in the superconducting 
state versus a for the model defined by Eqs.((TJ, (|58p . and 
(|59p with [7 = 0. Inset: the enlarged plot for small a 



Here we have ignored the corrections to the three-point 
vertex functions due to the superconducting gap again 
assuming that there is no strong ferromagnetic fluctua- 
tion. We are concerned with the case that the typical 
size of the spin-orbit splitting 2q;|<(p)| is much larger 
than |Afc|; i.e. the situation relevant to CePtaSi and 
CeRhSia. In this case, the second term of the right-hand 
side of piip . which is the van-Vleck-like contribution, 
is approximated by the expression for the normal state 
Xiix |53|) , and not affected by the superconducting tran- 
sition, while the Pauli term, i.e. the first term of (|llip . 
vanishes at zero temperature. When the energy depen- 
dence of the density of states is sufficiently small, e.g. 
the case of the spherical Fermi surface e = fc^/(2m), we 
have xl^(T = 0)/x^^(T = T^) = 1/2, as pointed out 
by Edelstein, and Gorkov and Rashba?^ We would like 
to stress that the nonzero x at zero temperature in the 
superconducting state is not due to the admixture of the 
spin singlet and triplet pairs, but rather attributed to 
the existence of the van-Vleck-like term caused by the 
spin-orbit interaction. 
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FIG. 10: Diagrams for the self-energy on the order of . 
Up to this order, the diagram (a) gives the most prominent 
electron correlation effects. The contributions from the other 
diagrams are suppressed by factors of order ^ 0{a\tQk\/ Ef) 
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FIG. 11: x^^{T = 0)/xxx{T = Tc) in the superconducting 
state versus U/W with a — 0.5 (open circles) and a — 1.0 
(crosses) calculated by the second order perturbative expan- 
sion with respect to U for the model defined by Eqs.([T|, (|58|l . 
and ((59)) . 



C. Implication for the Knight shift measurement of 
CePtsSi 

Here we would like to argue an implication of (|111[) 
for CePtaSi. There are strong pieces of experimental 
evidence suggesting the realization of an unconventional 
pairing state in CePtaSii^'SSii^'^ Nevertheless, the sym- 
metry of the Cooper pair has not yet been clarified so 
far. According to the Knight shift measurements done 
by Yogi et al^^^ for this system, both Xxx and Xzz does 
not show any significant change even below the supercon- 



ducting transition temperature Tc- The absence of the 
substantial decrease of Xxx for T < Tc seems not to be 
consistent with Eg. ljllip . which indicates the decrease of 
the Pauli term both in the spin singlet and triplet super- 
conducting states. To resolve this point, we note that the 
magnitude of the ratio Xxx{T — 0)/xxx{T — Tc) in the 
superconducting state crucially depends on the details of 
the electronic structure. For instance, when the density 
of states has the strong energy-dependence as shown in 
FIG. 7, the value oixxx{T = 0)/xxx{T = Tc) may deviate 
from 1/2. To demonstrate this, we use the model con- 
sidered in Sec. II. C, given by the Hamiltonian ^ with 
Eqs.linil) and ((59]) at the half-filhng. This model has the 
density of states similar to that depicted in FIG. 7 because 
of the existence of the van Hove singularity. We chose the 
magnitude of the superconducting gap as A± « 0.01, 
and calculate the spin susceptibility at zero tempera- 
ture from Eq. piip . In FIG. 8, we show the plot of 
XxxiT — 0)/XxxiT — Tc) versus the strength of the spin- 
orbit interaction a for the non-interacting case U — 0. 
It is seen that, for small a, Xxx{T = 0)/xxx{T = Tc) 
is nearly equal to zero. As a increases, the value of the 
ratio becomes larger, and eventually goes over 1/2, be- 
cause of the enhancement of the density of states due to 
the van Hove singularity. The maximum value of the ra- 
tio is still smaller than the experimentally observed value 
for CePtsSi, i.e. XxxiT = 0)/XxxiT = Tc) « 1. How- 
ever, it is expected that if electron correlation effects are 
included, the van-Vleck-like term enjoys the enormous 
enhancement caused by the contributions from the van 
Hove singularity, which is absent for the Pauli term if the 
spin-orbit splitted Fermi surfaces are sufficiently away 
from the van Hove points. To confirm this prediction, 
we carry out the perturbative calculation of the electron 
correlation effects by expanding the self-energy in terms 
of the Coulomb interaction U up to the second order. We 
consider the case where the spin-orbit splitting is much 
larger than the superconducting gap. Then, Xxx{T = 0) 
can be well approximated by the van-Vleck-like suscep- 
tibility in the normal state. Using the expressions in 
the normal state Eqs. (|52p and (j53p . we compute numer- 
ically the ratio XxxiT = 0)/Xxx{T = Tc) for a = 0.5 
and a = 1.0. Since the spin-orbit splitting is still much 
smaller than the band width w 8, we neglect cor- 
rections of order 0(Q;|tofe|/W^)j preserving only the con- 
tribution from the diagram FIG. 9(a). This amounts to 
considering only the term proportional to dT^x/ d{iJ.^Hx) 
in Eqs. (|55p and (j56p . Since the numerical computation 
of the derivative of the self-energy with respect to a mag- 
netic field is somewhat difficult, we use the approxima- 
tion d'Sx/ddiBHx) ~ d'S/d{ie), which is justified when 
a particular mode of spin fluctuation is not developed, as 
in the case of CePtaSii^ The calculated results of the ra- 
tio plotted against U/W are shown in FIG. 10. Here W is 
defined by the average of the band widths of ek+ and Sk- ■ 
As U/W increases, the enhancement of the van-Vleck-like 
susceptibility due to electron correlation overcomes that 
of the Pauli term, and the ratio Xxx{T — 0)/xxx{T = Tc) 
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approaches 1.0. Although the strength of the electron- 
electron interaction for which the ratio is nearly equal 
to 1.0 is considerably large, we believe that the result 
obtained by the perturbative calculation may be qualita- 
tively correct, and will be improved by including higher 
order corrections as long as the system is in the Fermi 
liquid state. Actually, CePtsSi is a Fermi liquid with 
the mass enhancement factor of order ~ 60, in which the 
spin susceptibility is also substantially enhanced. Since 
we cannot derive such a large mass enhancement from ab 
initio calculations, we regard the parameter U as an ef- 
fective interaction renormalized by correlation effects in 
the above calculations. 

Thus, it is possible that for CePtsSi, the absence of the 
decrease of Xxx (T) below Tc observed by the NMR exper- 
iment may be attributed to a specific energy dependence 
of the density of states, which leads the existence of the 
van-Vleck-like term much larger than the Pauli term. Un- 
fortunately, the LDA calculations of the electronic struc- 
ture of this system done by some groups are restricted 
within the case without the antiferromagnetic phase tran- 
sition, which occurs at T/v ~ 2K in CePtaSi.'^^'SS, Because 
of this reason, it is difficult to clarify from first principle 
calculations whether such a strong energy dependence of 
the density of states exists or or not for CePtaSi. Exper- 
imental studies which directly probe the energy depen- 
dence of the density of states for this system are desirable. 



D. Magnetoelectric effect: Paramagnetic 
supercurrent 

In the absence of inversion symmetry, an intriguing 
magnetoelectric effect exists in the superconducting state 
as well as in the normal state, as extensively argued by 
Edelstein, Yip, and the present author.— i^^'^'^i^^ In the 
case of the Rashba interaction with a potentail gradient 
along n = (0,0, 1), an in-plane magnetic field B induces 
a paramagnetic supercurrent, J = K,{n x B). Here K. 
is the magnetoelecric effect coefficient, of which the ex- 
plicit expression including electron correlation effects is 
presented in refiii. The inverse effect is also possible; 
i.e. the diamagnetic supercurrent flow gives rise to a 
magnetizationAi In the case of the Dresselhaus interac- 
tion with £o(fc) = {kx{kl - fc^), kyikl - k1),kz{kl - fc^)), 
a similar effect exists, but the paramagnetic supercurrent 
flows parallel to the applied magnetic field, J = ICB, in 
contrast to the Rashba case where the current flows per- 
pendicular to the magnetic field. The magnetoelectric 
coefficient K. is also directly related to the stability of the 
helical vortex phase considered by Kaur et al., and de- 
termines the magnitude of the center of mass momentum 
for inhomogeneous Cooper pairs."'^^ 

Recently, it was pointed out by Yip that the paramag- 
netic supercurrent is exactly canceled with the magneti- 
zation current Jm = cV x M in the complete Meissner 
state in a semi-infinite systemji^ However, in finite sys- 
tems, this cancellation is imperfect, and thus the param- 



agnetic supercurrent gives nonzero contributions to the 
bulk current. Moreover, in the mixed state, the cancella- 
tion is not complete generally. Here we would like to re- 
examine electron correlation effects on the paramagnetic 
supercurrent taking into account the partial cancellation 
with the magnetization current in the case of the Rashba 
interaction. 

Following Yip, we start from the following relations 
for the supercurrent density and the magnetization 
density M,W- 

= :A-(W0- — A) +/C(n X B), (112) 
2eA c 

^[nx (W0--A)]+Mzcc. (113) 
2e c 

Here the first term of the right-hand side of Eq. (|112p is 
the diamagnetic supercurrent density which is denoted 
by J'^''' = (Jf^, J^'^, Jf'^) in the following. The sec- 
ond term of the right-hand side of pi2p is the paramag- 
netic supercurrent induced by the inversion-symmetry- 
breaking spin-orbit interaction. The coefficient 1/A of 
the diamagnetic supercurrent is related to the super- 
fluid density, and equivalent to the Drude weight at zero 
temperature. The first term of the right-hand side of 
Eq. ()113p is the magnetization induced by the magneto- 
electric effect, and the second term Mzcc is a magne- 
tization caused by the usual Zeeman effect. For small 
B, Mzco = X-B- Eqs. pT2)) and (fTT3)) are obtained from 
the free energy in the absence of inversion symmetry as 
argued by several authorsi ^^'^'^i"'^^'^^ 

In addition to the supercurrent Jg, there exists the 
magnetization current density J m = cV x Af , which is 
derived by adding the term 

- J drM B = J dr{V x M) ■ A (114) 

to the free energy. Using Eqs. (|112p . (|113p . and the rela- 
tion V X J'^'^ = -B/(cA), we have^ 

K.{n X B) = -cV X (Af - Mzcc) 

+c/CA(-a, Jf ^ dx-Jf" + dyjf''). (115) 

Then, the total current density J^tot is recast into 

•^tot = J s + J M 

= J'^'"' + cV X Mzcc 

+ C}CK{^dxjf\ dyjf\ dxjf^ + dyjf^){l\&) 

The second term of the right-hand side of Eq. (|116p is the 
magnetization current caused by the usual Zeeman effect. 
The last term of Eq. (|116p is the paramagnetic supercur- 
rent due to the magnetoelectric effect with which we are 
concerned. Apparently, this paramagnetic term vanishes 
in the thermodynamic limit in the complete Meissner 
state in an infinite x-y planar system, as was pointed 
out by Yip." Nevertheless, it still plays a crucial role in 
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finite systems in the Meissner state, or generically in the 
mixed state. Using Ea. (|116p . we can discuss electron cor- 
relation effects on the paramagnetic supercurrent partly 
canceled with the magnetization current. According to 
the analysis presented in refJ^i^, /C is given by 
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(117) 



which is enhanced by the factor I/z^t, when the Wilson 
ratio is nearly equal to 2, as in the case of typical heavy 
fermion systems. On the other hand, the coefficient of the 
diamagnetic supercurrent 1 /A is suppressed by the factor 
ZkT ■ This fact indicates that the last term of the right- 
hand side of (|116p is enhanced by the factor 1/z^^ com- 
pared to the first term J*^^^. Thus, even when the partial 
cancellation with the magnetization current takes place, 
the paramagnetic supercurrent induced by the magneto- 
electric effect is still amplified in strongly correlated elec- 
tron systems with large mass enhancement. This obser- 
vation is quite important for the experimental detection 
of the magnetoelectric effect in heavy fermion supercon- 
ductors such as CePtaSi, CeRhSia, Ulr, and CelrSia. 



IV. SUMMARY 

In this paper, we have explored magnetic properties 
and transport phenomena in interacting electron systems 
without inversion symmetry in the normal state as well 
as in the superconducting state on the basis of the gen- 
eral framework of the Fermi liquid theory, taking into 
account electron correlation effects in a formally exact 
way. We have obtained the formulae for the transport 
coefficients related to the anomalous Hall effect, the ther- 
mal anomalous Hall effect, the spin Hall effect, and the 
magnetoelectric effect, which are the results of the par- 
ity violation. It is found that the spin Hall conductivity 
is not renormalized by electron correlation effects, and 
is determined solely by the band structure. Also, it is 



pointed out that in the superconducting state, the ther- 
mal anomalous anomalous Hall conductivity divided by 
the temperature and the magnetic field k^^^ /{TH^) re- 
mains finite even in the zero temperature limit, despite 
the absence of Bogoliubov quasiparticles carrying heat 
currents. This intriguing property is due to the fact that 
the anomalous Hall conductivity in the absence of in- 
version symmetry is dominated by electrons occupying 
the momentum space sandwiched between the spin-orbit 
splitted two Fermi surfaces, which are not affected by 
the superconducting transition. We have also examined 
electron correlation effects on the paramagnetic supercur- 
rent caused by the magnetoelectric effect taking account 
of the partial cancellation with the magnetization cur- 
rent. The experimental detection of these effects in the 
recently discovered noncentrosymmetric heavy fermion 
superconductors CePtaSi, Ulr, CeRhSia, and CelrSia. is 
an important future issue. 

Furthermore, we have demonstrated that in the normal 
state the temperature dependence of the van-Vleck-like 
spin susceptibility behaves like the Pauli susceptibility, in 
contrast to the usual van Vleck orbital susceptibility, and 
that the ratio of the van-Vleck-like term to the Pauli term 
depends crucially on the details of the electronic struc- 
ture and electron correlation effects. This result leads us 
a possible explanation for the recent NMR experimental 
data of CePtaSi, which indicates no change of the Knight 
shift below Tc for any directions of a magnetic field, con- 
trary to previous theoretical predictions.-® It is suggested 
that the strong energy dependence of the density of states 
and electron correlation effects may strongly enhance the 
van-Vleck-like susceptibility compared to the Pauli term, 
and thus the total spin susceptibility is not much affected 
by the superconducting transition. 
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